Degenerate Perturbation Theory Describing the Mixing of Orbital Angular 
Momentum Modes in Fabry-Perot Cavity Resonators 
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We present an analytic perturbation theory which extends the paraxial approximation for a com- 
mon cylindrically symmetric stable optical resonator and incorporates the differential, polarization- 
dependent reflectivity of a Bragg mirror. The degeneracy of Laguerre- Gauss modes with distinct 
orbital angular momentum (OAM) and polarization, but identical transverse order N, will become 
observably lifted at sufficiently small size and high finesse. The resulting paraxial eigenmodes possess 
two distinct OAM components, the fractional composition subtly depending on mirror structure. 
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Polarization-dependent effects in three-dimensional 
optical systems have in recent years received attention 
under the aspect of orbital angular momentum (OAM) 
It has become important to understand the OAM 
interactions with interfaces @, Q , waveguides Q and res- 
onators. One particular mechanism that warrants inves- 
tigation is that of small corrections to the paraxial theory 
of resonators giving rise to what may be regarded as op- 
tical spin-orbit coupling [f| . 

In this Rapid Communication, we extend electromag- 
netic resonator theory [ll| to provide a complete pertur- 
bation analysis and numerical computations for an opti- 
cal cavity which has an axis of rotational symmetry (z), 
but nevertheless does not conserve the component £ of 
OAM along that axis. This phenomenon itself is remark- 
able because the model system, shown in Fig. 1(a), ap- 
proaches the paraxial limit in which OAM conservation 
might be taken for granted. Standard paraxial modes 
may be chosen to have well-defined £ because the po- 
larization is transverse to z and factors out of the wave 
problem, leaving a scalar Helmholtz equation which maps 
to a quantum harmonic oscillator Q. Labeling the re- 
sulting transverse spectrum by £ € Z and a radial node 
number p 6 No, all modes with the same transverse or- 
der N = 2p + \£\ are degenerate in this set of approxi- 
mations. Therefore, going beyond this theory entails a 
degenerate perturbation theory for which we construct a 
coupling Hamiltonian V with the nominal Gaussian di- 
vergence angle, #d = kwo/2, as its small parameter. Here 
k is the wave number and wq is the waist radius. Since 
V is typically not diagonal in £, the solutions may be far 
from £-eigenstates. This holds even for arbitrarily small 
#d ■ Both polarization and £ mix in a paraxial spin-orbit 
coupling [g, Q. Our work, in part, yields a new way 
of generating OAM and other non-uniformly polarized 
light. Furthermore, in work involving cavity quantum 
electrodynamics, paraxial theory at the level developed 
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Figure 1: (a) Plano-concave model geometry. The top mirror 
radius of curvature R is 100 ^im in the examples, (b) Wave 
number splitting Afc = k<z — fci = Shi — Ski from numerical 
data at L = 25 /jm, A w 400 nm, indicating the avoided cross- 
ing between a mixable pair of modes, (c) Numerical results 
(symbols) for the mode mixing angle a closely fit the (solid) 
curve a — arctan[(e p — e s )/Ae — b]/2. Different symbols (col- 
ors) represent data series at L = 2.5 to 25 fim (in steps of 
2.5 /im) with kR » 780 to 1631. Ae and b are fit here for each 
series. 



here may be needed to distinguish spectral anti-crossings 
arising from passive cavity physics from those of strong 
photon-electron coupling. 

The main result of our work is an analytical, quantita- 
tive expression for the degree of OAM mixing in a dome- 
shaped cavity with a Bragg mirror. The latter is crucial 
to achieving high finesse in optical microcavities. Our 
results also shed new light on the recently-observed cou- 
pling between spatial and polarization degrees of freedom 
in broad-area vertical-cavity surface emitting lasers Q: 
polarization mixing was found to cause surprising spec- 
tral complexity even with simple (square) boundaries, ac- 
companied by intricate polarization patterns in the far- 
field emission. The essential physics is provided by the in- 
teraction of differently polarized plane-wave components 
with the planar, high-reflectivity distributed Bragg re- 
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fleeter (DBR) on which both Ref. [| and our system in 
Fig. 1 (a) are based. 

Typical DBRs comprise dielectric multilayers and have 
a "form birefringence": the reflectivities for TE (s- 
polarized) and TM (p-polarized) plane waves, r s ( P ) = 
|r s ( p ) | exp(«0 s ( p )), are unequal in phase at nonzero an- 
gles of incidence, 9. Previously [a, [Kj, we had found 
numerically that the reflection phase difference <p s (0u) — 
<Ap(0d) 7^ is responsible for the polarization mixing. 
A consequence of great practical importance is that by 
tuning cavity or mirror parameters, the coupling of OAM 
imposed by Maxwell's equations can be rigorously turned 
off for selected modes, giving them well-defined I. Aside 
from the consideration of OAM, our perturbation theory 
is fundamental to paraxial theory itself: we perform a 
significant extension to and completion of previous work 
by Yu and Luk [If] which derived the lowest order correc- 
tions to the paraxial modes of two-mirror cavities having 
perfect electrically conducting mirrors, r s = r p = — 1. 
Our approach, which allows a more general, dielectric 
planar mirror, requires the construction of a 2 x 2 per- 
turbation matrix V and predicts very different results 
when r s and r p are different functions. 

The work of Ref. [ll[ and its precursors (e.g. Ref. fl2T ]) 
appears to be forgotten in the current literature. This can 
be attributed to the difficulty of experimentally observ- 
ing the small spectral splittings caused by slightly non- 
paraxial perturbations. Agreement with a degenerate 
perturbation theory for r s ( p ) = — f was experimentally 
demonstrated using microwaves [f2j where, due to com- 
parable resonator size and wavelength, the mode spacing 
is large enough to resolve the lifting of the N + f-fold de- 
generacies labeled by N. However, with recent progress 
in miniaturization |8| , comparable size parameters are be- 
coming accessible to optical cavities similar to Fig. f (a). 

To write the splittings described in Ref. [Tfj in a form 
we can use for the following discussion, let us first re- 
view the unperturbed vector-field basis spanning the N- 
th transverse multiplet. These are the Laguerre- Gauss 
(LG) modes, written in polar coordinates as a product 
Mjy(p,(j>,z)& s of a scalar part Mf^, having orbital an- 
gular momentum £, and a circular polarization vector 
<r s = (x + siy) /a/2 , where s = ±1 is the spin degree of 
freedom. Expressions for the LG modes and their Bessel 
wave decompositions are given in Ref. [1 3j . Perturbations 
that preserve rotational symmetry around z will only cou- 
ple those M^(p, <fi, z)& s for which the total angular mo- 
mentum m = 1 + s around this axis is the same [Lsj . We 
henceforth consider N and to to be fixed, nonnegative 
parameters; we will briefly discuss later the allowance of 
to < 0. The non-conservation of OAM emerges here be- 
cause s can generally [1 7fl take two values, corresponding 
to the basis states M^f~ &\ and M]^ +1 <t_i. No symme- 
try prevents these states from coupling, unless m = in 
which case time-reversal invariance applies. 

Nevertheless, no such OAM coupling is found in 
Ref. [ffj , where perturbations merely split the wave num- 



bers of the pair {M^ l=Fl (T±i} by an amount 

Afc LG = k(M%+ l cr^) - fc(M--^0 = (f ) 

The calculation assumes a dome of vertical length L, top 
mirror radius of curvature i?, and JV P ) = — 1. We will 
now show that a more realistic model for r s ( p ) leads to ac- 
tual resonator modes ^N,m,i an d *&N,m,2 having a given 
to but forming a rotation of the LG basis pair by a mixing 
angle a £ (— ir/4, 7r/4): 

f^N, m ,i{a)\ _ /cos a -sincA / M 7 ^ 1 a 1 \ , , 
\*N, m A a )j V sin « c°sa J \M™ +1 a-i) ' [) 

The states at a = ±7r/4 are hybrid modes, one of which is 
predominantly (though not completely) composed of TM 
plane waves, the other being predominantly TE (T3. 
This approximate polarization separation of the hybrid 
modes yields an intuitive picture of mode mixing; the 
separation acts as a lever arm by which the Bragg mirror, 
having <p s ^ <p p , rotates the eigenmode basis away from 
the LG modes and toward the hybrid modes. 

The OAM-non-conservation described above is exhib- 
ited by the paraxial modes of two-mirror axisymmetric 
cavity resonators which f) are of sufficiently small size 
with respect to wavelength, 2) have sufficiently narrow 
resonance widths (low loss), and 3) have at least one 
mirror for which 4> s ^ <f> v (commercial dielectric mirrors 
meet this requirement). The first two requirements are 
essential in splitting the degeneracies of high order Gaus- 
sian modes, and Eq. ([T]) allows us to estimate whether 
this is possible. To perform the degenerate perturbation 
theory in the presence of property 3), we first develop 
the perturbation Hamiltonian, V, in the two-mode basis 
{M^ Tl cr±i}. V must be symmetric (V21 = V12), and for 
our purposes may be taken to be traceless (V22 = — Vn)- 
V then has eigenvectors v"i = ( f°lna ),V2 = ( cjSS) w ^ n 

a = (f/2)arctan(V 12 /l/ 22 ), (3) 

and eigenvalues Sku2) ~ -F(^i2 + ^S) 1 ^ 2 - ^ s system 
parameters are varied, the elements of V change and an 
anti-crossing of the hybrid modes emerges, cf. Fig. 1(b). 

All deviations from the paraxial limit must be con- 
sidered to lowest order in 9^, which for our cavity is 
6*2 = 2/[ky/L(R- L)]. The physical derivation of V is 
facilitated by noting that in the anti-crossing scenario, 
a = is equivalent to V12 = and hence corresponds 
to the assumptions underlying the known result Eq. |T]). 
Thus, Eq. (JTJ) should be reproduced by our model at 
a = 0. In wave number units, we therefore set 

Although Eq. ([T]) was derived for ideal-metal cavities, 
our more general DBR boundary conditions do not af- 
fect Eq. (U|), because they shift all LG modes with the 
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same N equally. To explain this, consider the penetra- 
tion depths <5£ s (p) of each of M^^a^i into the mir- 
ror layers. For a plane wave of s or p polarization with 
incident angle 9 at a DBR, 6L s r p ) — <f> s ( p )(k,6)/(2k). 
In order to capture the relevant material properties of 
the DBR, we neglect transmission and expand its reflec- 
tion phase in the plane wave angle of incidence, 9, as 
&(p)(M) « Mk) + e s ( P )(fc)6 |2 . 

Any vectorial mode \& can be decomposed into az- 
imuthally symmetrized plane waves (Bessel waves) by 
defining a "tilde" operator such that i& s ( p \6) essentially 
denotes the amplitudes of the TE(TM) plane waves of 
polar angle 9. The paraxial ty s ( p \9) is nonzero only near 
9 S3 0, with the averaged reflection phase of \& being 




\&(9)\ 2 U0) + \* P (8)\ 2 M0) Qd9, (5) 
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Figure 2: Mode-coupling width Ae from fit (symbols) com- 
pared to perturbative prediction (solid line), plotted as a func- 
tion of the cavity length L at fixed 7? = 100 ^m. Error bars 
were obtained from the fits of Fig. 1(c). Because the an- 
alytic equation 1|10[> is wavelength-independent, all data for 
A as 400 nm and A ~ 800 nm fall onto the same curve. 



with the normalization J(|* s | 2 + \^/ p \ 2 )9 d9 = 1. Spe- 
cializing to the LG modes with the abbreviation A±i = 
MJ2? Tl <x±i, circular polarization leads to |A! tl (6')| 2 — 
|A^. 1 (6')| 2 . The harmonic-oscillator nature of the trans- 
verse field entails that Eq. ((SJ), with the above expan- 
sion for </> s (p) , depends only on the mode order N but not 
on I. This carries over to the average penetration depth 
(^)a±i = (0)a±i /C^Os an d hence the effective cavity 
length L + (SL) A is identical for both modes A±i; this 
then implies equal spectral shifts, as claimed above. 

To obtain the off-diagonal element V12, we apply the 
same penetration-depth argument to the special case a = 
7r/4 where Eq. ([2} yields the hybrid modes. Their plane- 
wave amplitudes, x (9) and *^ p ^ 2 (0)> can be writ- 
ten purely in terms of the LG amplitudes A P tl (0) using 
Eq. ijSJ), and their splitting, Afc hybr id = k N ^ t2 - fc/v,m,i) 
is given by 



Afchybrid = ~(k/L) (6L N , m ,2 - $L N>mj x) (6) £) 

€p ~ €s f 9 3 A p +1 (9)A^ 1 (9)d9, (7) 



where the width of the crossover interval is given by 



where k is the unperturbed wave number. Therefore, 
we reach |Afch y brid| ^ IA^lgI if the form birefringence 
quantity, |e p — e s |, is made large. On the other hand, 
a — > 7r/4 implies | Vx2 j 3> | V221 5 so that the eigenvalues of 
V in this limit are Skin) ~ I fVi2- Setting 2Sk2 equal to 
Eq. JJJ) and performing the 9 integral, one obtains 



V12 = V21 



e p (k) - e s (fc) 



\) 2 -m?-f. (8) 



For some simple dielectric mirrors, e p — e s may be swept 
across zero by varying the cavity length across the nomi- 
nal L at which the mode pair of interest has unperturbed 
k equal to the design (center) wave number of the mirror, 
fed- The mixing angle a can then be written as 



a = (l/2)arctan{[e p (fc) - e B (fc)]/Ae}, 



(9) 



The formulas above complete the lowest order degener- 
ate perturbation theory for paraxial mode mixing. Inter- 
estingly, Ae is independent of wavelength: for fixed cavity 
geometry and mode labels N, to, modes of different lon- 
gitudinal node number (along z) will have different k but 
identical anti-crossing behavior when a is plotted versus 
e p — e s . This universal functional form provides a robust 
way of tailoring any desired mixing angle a. Most im- 
portantly under the aspect of OAM non-conservation, we 
can tune Eq. © to a — 0. 

Once the cavity linewidth for the relevant paraxial 
modes becomes less than the mode separation, Afc = 
2Sk 2 > m/(ikLR), the two modes ^N,m.j given by 
Eqs. ([U O mTj) are resolved at slightly different k (or 
L). Any excitation of the cavity would generally have 
non-zero overlap with these modes, and thus compli- 
cated mode patterns [f| [l3[ can be generated by simple 
excitation. The magnitude of a is zeroth order in 9^, 
and excursions near the asymptotic values can be seen 
in Fig. 1(c). The magnitude of the relative frequency 
splitting, however, is 0(#d) as #d — * 0. We note that 
observed modes will not be restricted to to > 0. Ax- 
ial symmetry creates an exact two-fold degeneracy of the 
vectorial LG basis modes under the transformation in 
which m, £, and s switch sign [l(|. The presence of the 
exact degeneracy [l8| does not "wash out" the generation 
of complicated mode patterns and is more fully discussed 
in Refs. 0, 0,(11. 

We have numerically calculated modes for R = 100 /im 
and L — 2.5 to 25 /xm with a DBR comprising 36 pairs 
of quarter- wave dielectric layers A and B with refractive 
indices tia = 3.52 and Ub = 3.00, with layer B at the 
top surface. Typical values of e s ( p ) were around —3, with 
d(e p — e s )/d(fc — fed) ~ — 2.8 zxm. Data were analyzed 
for modes close to Ad = 2-7r/fcd « 400 and 800 nm. We 
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considered mode pairs with N = 2 and m = 1, the lowest 
values for which OAM mixing can occur. 

The numerical data for the mixing angle are fit ex- 
tremely well by Eq. §§§ if we allow for an offset b in the 
argument of the arctan, as done in Fig. 1(c). The com- 
parison between numerical fit and Eq. (I10|l is shown in 
Fig. 2. The offset b has median —0.32 for our data and 
empirically behaves as C/[kL(l — L/R)], where C is a 
slowly- varying function of k and fed- Taking the limit 
$rj — > such that L/R is bounded away from and 1 
implies that \b\ = 0(9^). This ncxt-highest-order cor- 
rection to our perturbation theory will be discussed in a 
subsequent publication [1 51 ] - 

To spectrally resolve the transverse mode splitting 
along the entire mixing curve, the cavity must obey 
2fci?[l— (R1R2) 1 / 2 ] < to, where Ri(2) are the power reflec- 
tivities of the two mirrors. Microwave experiments may 
be the most direct approach. Alternatively, paraxial spin- 
orbit coupling may be realized in microcavity resonators 
for quantum- information applications, cf. Ref. Q, where 
small size and high finesse are required. Such cavities 
could generate OAM or hybrid beams at light levels from 
single-photon-on-demand to that of a macroscopic laser. 
In particular, one could utilize the fine structure and var- 
ied spatial patterns of the split modes. This has partic- 
ular potential for quantum information applications: the 



order TV family provides N + 1 nearly degenerate energy 
levels corresponding to modes having different vectorial 
spatial patterns of the electric field. This spectral and 
spatial structure combined with quantum dots at the pla- 
nar mirror may possess quantum logic capability. 

In conclusion, both the mixing angle and the frequency 
splitting for OAM-mixed paraxial resonator modes in 
an axisymmetric cavity have been analytically derived 
here in a degenerate perturbation theory which includes 
the form birefringence of a practical mirror. The phe- 
nomenon is similar in principle to the polarization cou- 
pling observed in Refs. [9|, [l6| : differences in the penetra- 
tion depths for TE and TM plane waves modify the vec- 
torial resonator modes. For our case however, eigenmode 
coupling persists from the non-paraxial regime to the 
deeply paraxial regime. By examining the latter, we have 
shown here that the coupling can in fact be turned on and 
off via the material parameters e s ( p ). The underlying in- 
adequacy of a scalar paraxial treatment is washed out in 
macroscopic cavities, but must be regarded as a funda- 
mental limitation in high-finesse microcavities, where the 
correct starting point for any paraxial formulation must 
be the mode basis of Eq. which is heterogeneous in 
orbital angular momentum. 

This work is supported in part by National Science 
Foundation Grant No. ECS-0239332. 
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